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Measurements with an atomic force microscope (AFM) offer a direct way to probe elastic prop- 
erties of lipid bilayer membranes locally: provided the underlying stress-strain relation is known, 
material parameters such as surface tension or bending rigidity may be deduced. In a recent ex- 
periment a pore-spanning membrane was poked with an AFM tip, yielding a linear behavior of the 
force-indentation curves. A theoretical model for this case is presented here which describes these 
curves in the framework of Helfrich theory. The linear behavior of the measurements is reproduced 
if one neglects the influence of adhesion between tip and membrane. Including it via an adhesion 
balance changes the situation significantly: force-distance curves cease to be linear, hysteresis and 
nonzero detachment forces can show up. The characteristics of this rich scenario are discussed in 
detail in this article. 

PACS numbers: 87.16.Dg,68.37.Ps,02.30.Hq 



I. INTRODUCTION 

Lipid bilayer membranes constitute one of the most 
fundamental components of all living cells. Apart from 
their obvious structural role in organizing distinct bio- 
chemical compartments, their contributions to essential 
functions such as protein organization, sorting, or sig- 
nalling arc now well documented In fact, their tasks 
significantly exceed mere passive separation or solubiliza- 
tion of proteins, since often mechanical membrane prop- 
erties are intricately linked to these biological functions, 
most visibly in all cases which go along with distinct 
membrane deformations, such as exo- and endocytosis 
0, vesiculation 1, 0, @, viral budding 0, cytoskele- 
ton interaction [7], and cytokinesis Consequently, a 
quantitative knowledge of the material parameters which 
characterize a membrane's elastic response - most no- 
tably the bending modulus k — is also biologically desir- 
able. 

Several methods for the experimental determination 
of k have been proposed, such as monitoring the spec- 
trum of thermal undulations via light microscopy [9l.ll0||. 
analyzing the relative area change of vesicles under mi- 
cropipette aspiration [III, [l2j , or measuring the force re- 
quired to pull thin membrane tethers [H, 03, HE Ell- 
With the possible exception of the tether experiments, 
these techniques are global in nature, i. e., they supply 
information averaged over millions of lipids, if not over 
entire vesicles or cells. Yet, in a biological context this 
may be insufficient For instance, membrane proper- 
ties such as their lipid composition or bilayer phase (and 
thus mechanical rigidity) have be en p roposed to vary on 
submicroscopic length scales [3, Il9j . Despite being bi- 
ologically enticing, this suggestion, known as the "raft 
hyp othesis" , has repeatedly come under critical scrutiny 
[20l . I2H , HU , precisely because the existence of such small 



domains is extremely hard to prove. 

An obvious tool to obtain mechanical information for 
small samples is the atomic force microscope (AFM) 
[HI, and it has indeed been used to probe cell elastic 
propert ies (such as for instance their Young modulus) 
[2J,|25(. Yet, obtaining truly local information still poses 
a formidable challenge. Apart from several complications 
associated with the inhomogeneous cell surface and intra- 
cellular structures beneath the lipid bilayer, one partic- 
ularly notable difficulty is that the basically unknown 
boundary conditions of the cell membrane away from the 
spot where the AFM tip indents it preclude a quanti- 
tative interpretation of the measured force, i. e. a clean 
way to translate this force into (local) material proper- 
ties. To overcome this problem Steltenkamp et al. have 
recently suggested to spread the cell membrane over an 
adhesive substrate which features circular pores of well- 
defined radius [26| . Poking the resulting "nanodrums" 
would then constitute an elasto-mechanical experiment 
with precisely defined geometry. Using simple model 
membranes, the authors could in fact show that a quan- 
titative description of such measurements is possible us- 
ing the standard continuum curvature-elastic membrane 
model due to Canham [27| and Helfrich (28j . 

Spreading a cellular membrane without erasing inter- 
esting local lipid structures obviously poses an experi- 
mental challenge; but the setup also faces another prob- 
lem which has its origin in an "elastic curiosity" : even 
significant indentations, which require the full nonlin- 
ear version of the Helfrich shape equations for their cor- 
rect description, end up displaying force-distance-curves 
which arc more or less linear — a finding in accord with 
the initial regime of membrane tether pulling (29l . [30| . 
Yet, this simple functional form makes a unique extrac- 
tion of the two main mechanical properties, tension and 
bending modulus, difficult. Is the nanodrum setup thus 
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FIG. 1: Illustration of the geometry. A parabolic tip with cur- 
vature radius J? t i P indents a pore-spanning membrane with a 
force F to a certain depth ho- The radius of the pore is 



The membrane detaches from the tip at a radial dis- 



tance p = c. The two possible parametrizations h(p) and i/>(s) 
are explained in the beginning of Chapter Mil 



of which depend in a predictable way on the experimental 
setup. In our case we will for instance find snap-on and 
snap-off events between tip and membrane, which rest on 
the fact that binding is not pre-determined, and whose 
correct description is very important for reliably inter- 
preting any AFM force experiments. Moreover, we will 
also see that the very occurrence of tethers is a much 
more subtle phenomenon, since an adhering membrane 
pulled upwards may in fact prefer to detach rather than 
being pulled into a tether — a question treated previously 
(and on the linear level) by Boulbitch (44[. 

Our paper is organized as follows: in Chapter [TT] we 
introduce the model of our system and discuss the rele- 
vant energies. In Chapter IIIII we present the equations 
that have to be solved in order to find membrane profiles, 
force- indentation curves and detachment forces. This will 
also include a treatment of the nonlinear case which was 
only mentioned very briefly in Ref. [26| . In Chapter IIVI 
the results of our calculations are summarized and com- 
pared to existing [26[ measurements. We end in Chap- 
ter [V] with a discussion how the predictions for indenta- 
tion and adhesion characteristics can be used to extract 
material properties in future experiments. 



futile? 

In the present work we develop the theoretical ba- 
sis for a slight extension of the nanodrum experiment 
that will help to overcome this impasse. We will 
show that an additional adhesion between the AFM 
tip and the pore-spanning membrane will change the 
situation very significantly — quantitatively and qualita- 
tively. Force-distance-curves cease to be linear, hystere- 
sis, nonzero detachment forces and membrane overhangs 
can show up, and various new stable and unstable equi- 
librium branches emerge. The magnitude and charac- 
teristics of all these new effects can be quantitatively 
predicted using well established techniques which have 
previously been used successfully to study vesicle shapes 
[HI H3, [HiJMJpl, vesicle adhesion j3pjpL colloidal 
wrapping (MSI!! or tether pulling [MM 5J 51 S3 • 



[53 or tether pulling [29|,|3C 
The key "ingredient" underlying most of the new 
physics is the fact that the membrane can choose its 
point of detachment from the AFM tip . Unlike in the 
existing point force descriptions [2^. |30||. in which a cer- 
tain (pushing or pulling) force is applied at one point of 
the membrane, our description accounts for the fact that 
the generally nonvanishing interaction energy per unit 
area between tip and membrane co-determines their con- 
tact area over which they transmit forces, and thereby 
influence the entire force-distance-curve. What may at 
first seem like a minor modification of boundary condi- 
tions quickly proves to open a very rich and partly also 
complicated scenario, whose various facets may subse- 
quently be used to extract information about the mem- 
brane. In fact, Smith et al. [HI, S3] have demonstrated in 
a related situation that the competition between adhe- 
sion and tether pulling for substrate-bound vesicles gives 
rise to various first- and second-order transitions, details 



II. THE MODEL 



A. Geometry of the system 



We consider a flat solid substrate with a circular pore 



of radius i? D 



A lipid bilayer membrane is adsorbed 



onto the substrate and spans the pore. In the situation 
we want to analyze an AFM tip is used to probe the prop- 
erties of the free pore-spanning membrane. We assume 
that the tip has a parabolic shape with curvature radius 
i?tip at its apex. Furthermore, we restrict ourselves to 
the static axisymmetric situation in which the tip pokes 
the free-standing membrane exactly in the middle of the 
pore (see Fig. [T]). 

For a certain downward force F > the membrane 
is indented to a corresponding depth ho > which is 
measured from the plane of the substrate to the depth 
of the apex of the tip. Note that it is also possible to 
pull the membrane up with a force F < in the opposite 
direction if attractive interactions attach the membrane 
to the tip. 

In the following, we will model the bilayer as a two- 
dimensional surface. This is a valid approach provided 
the thickness of the membrane (approx. 5 nm) is much 
smaller than (i) the membrane's lateral extension as well 
as (ii) length scales of interest such as local radii of cur- 
vature. 

With this geometric setup in mind, let us now consider 
the different energy contributions we want to include in 
our model. 
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B. Energy considerations 

The total energy of the system "pore-tip" comprises 
different contributions: the membrane is under a lateral 
tension a. To pull excess area into the pore, work has to 
be done a gain st the adhesion between membrane and flat 
substrate |45l|. It is given by a times the excess area [46j . 
Additionally, a curvature energy is associated with the 
membrane. According to Canham [27| and Helfrich [28| 
the Hamiltonian for an up-down symmetric membrane is 
then 



E, 



last 



dA (a 



-K z + RKc 



(1) 



where £ denotes the surface of the membrane part which 
spans the pore. The proportionality constants k and R 
are called bending rigidity and saddle-splay modulus, re- 
spectively. The Gaussian curvature Kq is the product 
of the two principal curvatures whereas K is their sum 
[47l 148} . Note that the last term of energy (JTJ) yields zero 
in our specific problem [49(. 

With the help of the two material constants a and k 
one can define a characteristic lengthscale 



(2) 



which does not depend on geometric boundary conditions 
such as the radius of the tip or the pore but only on 
properties of the membrane. On scales larger than A 
tension is the more important energy contributions; on 
smaller scales bending dominates. 

Apart from tension and bending, an adhesion between 
tip and membrane may contribute to the total energy. 
We assume that it is proportional to the contact area 
^contact between tip and membrane with a proportional- 
ity constant w, the adhesion energy per area. 

If the indentation ho is given and one wants to deter- 
mine the force F, the total energy can thus be written 
as 



-K 2 



W^4contact 



(3) 



Under certain circumstances, however, it is more conve- 
nient to consider the problem for a given force F. Both 
ensembles ("constant indentation" vs. "constant force") 
are connected via a Legendre transformation (4l|, [4^ |. 
^totai = ^totai — While the ground states one ob- 

tains for the two ensembles will be the same, questions of 
stability depend on the ensemble: a profile found to be 
stable under constant height conditions is not necessarily 
stable under constant force conditions. 

The route we want to follow here in order to find force- 
indentation curves is to determine the equilibrium shapes 
of the non-bound section of the membrane via a func- 
tional minimization. The energy contributions caused 
by the bounded section of the membrane enter via the 
appropriate boundary conditions (see Chapter IIIII and 



Appendix [5J ■ These imply that the contact point c is 
not known a priori but has to be determined as well 
("moving boundary problem"). In the next section we 
will show how one can set up the appropriate mathemat- 
ical formulation of the problem to get membrane profiles 
and force-indentation curves. 



III. SHAPE EQUATION AND APPROPRIATE 
BOUNDARY CONDITIONS 



To describe the shape of the membrane we use two dif- 
ferent kinds of parametrization (see Fig.[T]): for the linear 
approximation it is sufficient to use "Monge" gauge where 
the position of the membrane is given by a height h(p) 
above (or below) the underlying reference plane. The 
disadvantage of this parametrization is that it does not 
allow for "overhangs" . Since these may be present in the 
full nonlinear problem, we will use the "angle-arclength" 
parametrization in the exact calculations: the angle ip(s) 
with respect to the horizontal substrate as a function of 
arclength s fully describes the shape. 



A. Linear approximations 

To get the profile of the free membrane one has to 
solve the appropriate Euler-Langrange ("shape") equa- 
tion. This equation is typically a fourth order nonlinear 
partial differential equation and thus in most cases im- 
possible to solve analytically. One may, however, con- 
sider cases where the membrane is indented only a little 
and gradients are small. In that case one may linearize 
the energy functional. In the constant indentation en- 
semble one gets for the free part 



E = 



dAn 



(4) 



where d^4|| is the area element on the flat reference 
plane and Ef ree is the projected surface of the free pore- 
spanning membrane. The symbol V denotes the two- 
dimensional nabla operator in the reference plane. 

The appropriate shape equation can be derived by set- 
ting the first variation of energy ([4]) to zero, yielding 



V 2 (V 2 - \- 2 )h = 



(5) 



The solution to this equation is a linear combination of 
the eigenfunctions of the Laplacian corresponding to the 
eigenvalues and A -2 . For axial symmetry it is given by 
h(p) = hi + h 2 In (p/X) + h 3 I {p/X) + h i K (p/X), where 
Jo and K are the modified Bessel functions of the first 
and the second kind, respectively The constants 

hi, . . . , hi are determined from the appropriate boundary 
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conditions (see Appendix [5} : 



as [33|,134|, 135, 39, 41]: 



h(Rporc) 
h (Rpore) 

and h"(c) 



, h(c) = -h a + 
, h'(c) 



2R 



tip 



-Rtip 



R 



lip 



(6a) 
(6b) 

(6c) 



where a dash denotes a derivative with respect to p. 
Even though the differential equation is of fourth order, 
five conditions are required due to its moving bound- 
ary nature, i. e., c is to be determined from an adhesion 
balance — which is in fact the origin of the fifth condition 
(f6"c|) (see Appendix [AJ . 

The solution of the boundary value problem (|5l6p can 
be used in two ways to calculate the force for a prescribed 
indentation: first, one can insert the profile back into the 
functional ^ to obtain the energy of the equilibrium 
solution, which will then parametrically depend on the 
indentation ho- Its derivative with respect to ho yields 
the force F. Second, one can also consider stresses: in 
analogy to elasticity theory JM| F is given by the integral 
of the flux of surface stress [52| through a closed contour 
around the tip. 

The second approach is used in the present work; it has 
the advantage that the final expression for the force can 
be written in a closed form [26| (see also Appendix IB")) : 



F 



dK 



P=R P 



(7) 



This equation is exact. Inserting the solution h(p) of the 
boundary value problem (|5I6|) into (O yields the value of 
the force in the linear regime. 

A little warning might be due here: expression ([7]) is 
evaluated at the rim of the pore where the profile is 
flat even for high indentations. One might thus won- 
der whether inserting the small gradient solution would 
actually lead to an exact result. This is, however, not 
the case, because the membrane shape at the rim pre- 
dicted by the linear calculation is not identical with the 
prediction from the full nonlinear theory — except for its 
flatness, which is enforced by the boundary conditions. 
There is no magical way to avoid solving the nonlinear 
shape equation if one wants the exact answer. 



B. Complete nonlinear formulation 



E 



As L 



d.s 



sin-0 



2p 
A 2 



+ X p (p — cost/;) + \ z (z — svaip) — Fz\ , (8) 

where s is the arclength at the contact point c and s the 
arclength at -R p0 rc- The dot denotes the derivative with 
respect to s. The Langrange multiplier functions A p and 
X z ensure that the geometric conditions p = cos ip and 
z = sin?/; are fulfilled everywhere. 

In order to make the numerical integration easier let 
us rewrite the problem in a Hamiltonian formulation [33l . 
H3, [H, |4l|: the conjugate momenta are — dL/dtp = 
2p[ip + sm(ip)/p], p p = dL/dp — A p , and p z = dL/dz — 
X z — F. The (scaled) Hamiltonian is then given by 



H = ijjp,p + pp p + zp z - L 



_ Py> sin ip _ 2p 

~ P p a7 



+ p p cos ip + {p z + F) sin ip. (9) 



Note that H is not explicitly dependent on s and is thus 
a conserved quantity. Instead of one fourth order one 
then has six first order ordinary differential equations, 
the Hamilton equations: 



?± _ 
2p 

cosip 



— — = simp 



Pw 



Pp 



dH 
dp^ 
dH 
dp P 
dH 
dp 
dH 

9-0 
dH 
~~dp~ 
dH 

dz 



sini/) 



P 

p±fP± _ 

P 
. 



(10a) 
(10b) 
(10c) 

+ F) cos ip + p p sin tp (10d) 
sintp^ 



P 



2 
A 2 



(lOe) 



(lOf) 



According to the last equation, p z has to be constant 
along the profile. Its value can be found by considering 
the integral over the flux of surface stress which has to 
equal the applied force. This implies that p z vanishes 
everywhere (see Appendix [El . 

Equations (flC)]) can be solved numerically subject to 
the boundary conditions (see also Appendices lAl and [C|l : 



Let us now shift to the angle-arclength parametrization 
and consider the full nonlinear problem. In principle, the 
constant height ensemble could be used here as well. It is, 
however, technically much easier to fix F instead in order 
to reduce the number of boundary conditions one has to 
fulfill at the rim of the pore (see below and Appendix [C"jl . 

In this paragraph all variables with a tilde are scaled 
with 7tk, i. e.: E := E/(ttk), F :— F/(ttk), etc. The en- 
ergy functional of the free membrane can then be written 



and H 



, ip(s) = a , 
(cos a) 3 / 2w 

"V7 



R 



tip 



o 



(lla) 
(lib) 
(11c) 



where contact radius c and contact angle a are connected 
via c = i?tiptana. The solution to (fTTJl fTTj) gives the 
indentation /iq for some prescribed force F. 
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IV. RESULTS 

This chapter will summarize the characteristic features 
of the solution to the boundary value problems {5j [6]) and 
(fTTH ITT]) . In addition, the theory will be shown to be in 
accord with available experimental results. 

We will introduce some additional variable rescaling 
in order to make generalizations of the results easier: 
lengths will be scaled with i?ti P - We also define 



2wRl p 



and / 



FR 



tip 



In a typical experiment the curvature of the tip is of the 
order of ten nanometer (5-40 nm) and pore radii may 
lie between 30 and 200 nm [53|. The bending rigidity 
of a fluid membrane may vary between one and a hun- 
dred IcbT [H3]. One expects a maximum surface tension 
of the order of a few mN/m, which is approximately the 
rupture tension for a fluid phospholipid bilayer [551 ] . A 
maximum value of the adhesion can be found by assum- 
ing that a few k^T per lipid is stored if membrane and 
tip are in contact. One arrives at w max w 10mJ/m 2 . For 
the continuum theory to be valid Eqns. (|6c[ lllb[) imply 
that \/2w max / k < 1/d, where d ~ 5nm is the bilayer 
thickness. This estimate yields approx. the same max- 
imum value for w max as before since k is at most 100 
k B T. 

Thus, <7 and w can in principal vary between and 
10 3 . Realistically, if we set i? t i P = 10 nm and consider 
a typical fluid phospholipid bilayer with k ~ 20 k^T, cr 
and w are of the order of 1. Furthermore, we will focus 
on a pore radius of i? p0 rc = 3 in the following. 

In order to understand, how adhesion energy modifies 
the force-distance behavior, let us first briefly revisit the 
case where there is no adhesion between tip and mem- 
brane (w = 0). 



A. No adhesion between tip and membrane 

In Fig. [2] the shapes of the membrane for different val- 
ues of indentation are presented in scaled units. The 
linear calculations are dotted whereas the exact result is 
plotted with solid lines. For small indentations the two 
solutions overlap; for increasing ho, however, the devia- 
tions become larger just as one ex pec ts for a small gradi- 
ent approximation (see also Ref. [39( for another exam- 
ple). While the differences are noticeable, they appear 
fairly benign, such that one would maybe not expect big 
changes in the force-distance behavior. We will soon find 
out that these hopes will not be fulfilled. 

A deeper indentation also means that the tip has to 
exert a higher force. In Figs. [3] and UJ log-log plots 
of force-distance curves for different values of cr are 
shown. The dashed line marks the maximum indenta- 




FIG. 2: Membrane profiles for different indentations ho, all for 
cr — 1 (solid lines: nonlinear calulations, dashed lines: linear 
approximation, grey shades: AFM tips). The corresponding 
forces /(ho) for the three different indentations are (nonlinear 
calculations): /(0.3) = 0.81, /(l) = 2.45, /(2) = 4.27. 



of tip and pore. In the limit of high forces all curves con- 
verge and approach /io,max; for small forces the curves are 
linear in /. Let us quantify the indentation response by 
defining the (scaled) apparent spring constant JC of the 
nanodrum-AFM system via 



K 



dh 



(13) 



tion ho i 



i?p OI . /2 which is allowed by the geometry 



A linear force-distance-curve has a constant JC and thus 
follows an apparent Hookean behavior / = JCho- In 
unsealed units, the spring constant is given by JC = 
OF/ dho — JCnn/R^. For typical values i? t i P = 10 nm 

and k = 20/cbT this implies K = K x 2.6 mN/m. 

The smaller a, the less force has to be applied to reach 
the same indentation (see Fig. [3]). For decreasing a the 
force-distance curves converge to the limiting curve of 
the pure bending case, for which a = 0; this is plotted 
dashed-dotted in Fig. [31 In the opposite pure tension 
limit (k — > or a — > oo) the curves become essentially 
linear in a, as can be seen clearly after scaling out the ten- 
sion (see Fig. UJ. It is possible to calculate this second 
limiting curve in the linear regime: the linearized Eu- 
ler Lagrange equation reduces to the Laplace equation, 
Ah = 0, which is solved by h(p) — d\ + d 2 ln(p/i? poro ) 
in the present axial symmetry. The constants d\ and d 2 
can be determined with the help of the two boundary 
conditions h(R polc ) — and h(c) — — ho + c 2 /2R t i p . The 
contact point c is then determined by a straightforward 
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FIG. 3: Force-distance curves for w — and a = jg, j, 1, 4, 16 
and 64 (a increasing from left to right). The curve for a = 
is dashed-dotted. The inset shows the corresponding scaled 
apparent spring constant K, (see Eqn. (JT3J) in the small force 
limit, illustrating its two different regimes of small and large 
tension with a crossover around a ~ 1. 
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FIG. 5: Comparison between experiment (solid grey line) 
and theory (dashed line: linear approximation; solid black 
line: nonlinear calculations) [26|. The theoretical curves 
are obtained with the following parameters: -R p0 rc = 90 nm, 
Rtip = 20 nm, k = 10" 19 J « 24 k B T, a = 1.1 mN/m. 
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FIG. 4: Scaled force-distance curves for w = and a = 
|,1,4, 16 and 64 (a increasing from right to left). The 
solution for a — > oo in the linear regime is dashed-dotted. 
Nonlinear results are plotted with solid lines, the linear ap- 
proximation is dotted. 



energy minimization. The final result for the indentation 
depth is: 



"0 



1 - In 



(AIL. 

\ ixt pore 



(14) 



which is plotted dashed-dotted in Fig. [4] At any given 
penetration the force is now strictly proportional to the 
tension. Notice also the remarkably weak (logarithmic) 
dependency of penetration on pore size. 

All force-distance curves presented in this section ex- 
hibit a linear behavior for small forces. In this limit the 
scaled spring constant for the systems just discussed is 



well described by the empirical relation JC ~ 1.76 + fx 0,89 
(see inset in Fig. [3]). Combining this with our observation 
that for typical system parameters K, = K X 2.6 mN/m, 
we see that a nanodrum's stiffness can be very well 
matched by available (soft) AFM cantilevers, showing 
that the suggested experiments are indeed feasible. In 
fact, Fig. [5] shows the results of such an indentation ex- 
periment (solid grey line). Here, a fluid DOTAP (1,2- 
dioleoyl-3-trimethylammonium-propane chloride) mem- 
brane was suspended over a pore of radius i? por c = 90 nm 
and subsequently probed with a tip of radius i?tip = 
20 nm [26| . The apparent spring constant is found to be 
3.9 mN/m. To fit the data we optimized the material pa- 
rameters a and k. The linear approximation (asymptot- 
ically) matches the curve down to an indentation depth 
of about 40 nm as one can see in Fig. [5] (dashed line). For 
larger indentations the small gradient assumption breaks 
down. The nonlinear calculation (solid black line) de- 
scribes the data correctly down to a much deeper pen- 
etration depth of 150 nm but diverges for larger values. 
This deviation is most likely not a failure of the elastic 
model but a consequence of our simplified assumptions 
for the tip geometry. As shown in Fig. IB of the Supple- 
mentary Information to Ref. [26[ the tip is parabolic at 
its apex, but further up it narrows quicker and assumes 
a more cylindrical shape. It therefore can penetrate the 
pore much deeper than one would expect if the parabolic 
shape were correct for the entire tip. 

Apart from this difficulty, theory and experiment are 
in good agreement. There is, however, a catch. Since 
we cannot trust the force-distance behavior close to the 
depth-saturation (due to its displeasingly strong depen- 
dence on the actual tip shape), the remaining inter- 
pretable part of the force-distance-curve is linear, and 
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its slope is the only parameter that can be extracted 
from the data [5(|. For the theoretical calculation one 
needs two parameters, a and k. Fitting both to a line 
is not possible. In Ref. [2|| this obstacle was overcome 
by estimating n from other measurements to be about 
1CP 19 J. The surface tension a could then be adjusted to 
l.lmN/m to match the data — which, reassuringly, is a 
very meaningful value. 

Alternatively, one may proceed in a different manner. 
In the experiment a small snap-off peak could be observed 
upon retraction of the AFM tip which was due to the at- 
traction between tip and membrane. Although this could 
be neglected in the interpretation of the measurements 
of Ref. [1^ , one may think of deliberately increasing the 
adhesion between membrane and tip in a follow-up ex- 
periment by chemically functionalizing the tip. With this 
additional tuning parameter one may get further infor- 
mation on the values of the material parameters in ques- 
tion. 
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FIG. 6: Force-distance curves for a — 1 and w = 0,1, 2, 3,4, 5 
(from right to left). The region of hysteresis in the curve for 
w = 2 is magnified in the inset. In this case the energy barrier 
at / = 0.414 is approx. 1 n. Overhang branches are omitted. 



B. Including adhesion between tip and membrane 

In the following, we will also allow for adhesion be- 
tween tip and membrane, i. e. w is not necessarily equal 
to zero. This will change the qualitative behavior of the 
force-distance curves dramatically: for fixed a and w dif- 
ferent solution branches can be found. A hysteresis may 
occur as well, as we will see in the next section. Addi- 
tionally, stable membrane profiles exist even if the tip is 
pulled upwards. It is therefore possible to calculate the 
maximum pulling force that can be applied before the tip 
detaches from the membrane and relate it to the value of 
the adhesion between tip and membrane. 



1. Weak adhesion energy 

In this section, we will first investigate the case of 
weak adhesion, w < 5. The scaled surface tension a 
will be fixed to 1. It turns out that once the tip is adhe- 
sive, "overhang" profiles may occur, i. e., shapes where 
at some point \ifj(s)\ > 90°. We will first ignore these 
solution branches and come back to them later. 

Fig. [6] illustrates force-distance-curves for ii = 
0, 1, . . . , 5. Compared to the nonadhesive case, for which 
an essentially linear behavior levels off towards maximum 
penetration, adhesive tips behave quite differently. Al- 
ready for w = 1 an initial Hookean response at small 
forces is soon followed by a regime in which the system 
displays a much greater sensitivity towards an externally 
applied stress, i. e., where the scaled spring constant K, 
drops at intermediate penetrations. Physically this of 
course originates from the fact that adhesion helps to 
achieve higher penetrations, because the tip is pulled to- 
wards the membrane, but notice that this does not lead 
to a uniform reduction of /C: softening only sets in be- 
yond a certain indentation. 



Shortly beyond w = 1 a point is reached where the 
force-distance-curve displays a vertical slope at which the 
apparent spring constant K, vanishes. For even larger 
values of adhesion a hysteresis loop opens, featuring a 
locally unstable region with JC < 0. This is the case for 
w = 2, and the region around the instability is magnified 
in the inset of Fig. [Sj Notice that the dotted branch cor- 
responding to JC < still belongs to solutions for which 
the functional ^ is stationary, yet the energy plotted 
against penetration ho (or, alternatively, contact angle 
a) has a local maximum, confirming that these solutions 
are unstable against contact point variations. The two 
dashed branches in the inset of Fig. [5] have a positive 
JC and correspond to local minima in the energy, how- 
ever, they are globally unstable against the alternative 
minimum of larger or smaller ho. The true global mini- 
mum is indicated by the bold solid curve, which exhibits 
a discontinuity at / = 0.414. 

Depending on the current scanning direction this hys- 
teretic force-distance-curve manifests itself in a snap-on 
or snap-off event. Such a behavior is reminiscent of a 
buckling transition (such as for instance Euler buckling 
of a rod under compression [5l| ) — with two caveats: first, 
notice that the membrane does not stay flat up to a crit- 
ical buckling force at which it suddenly yields; rather, 
the system starts off with a linear stress-strain relation 
and only later undergoes an adhesion-driven discontinu- 
ity. Appreciating this point is quite important for the 
interpretation of measured force-distance curves: upon 
approach of tip and membrane the snap-on will occur 
neither at zero force nor at zero penetration. Second, 
one should not forget that hysteresis is ultimately a con- 
sequence of the energy barrier which goes along with 
such discontinuities. For macroscopic systems this bar- 
rier is typically so big that the transition actually hap- 
pens at either of the two end-points of the S-shaped hys- 
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teresis curve, where the barrier vanishes (the "spinodal 
points"). However, for nano-systems barriers are much 
smaller, comparable to thermal energy k^T, such that 
thermal fluctuations can assist the barrier-crossing event. 
In the present case the barrier at the equilibrium tran- 
sition point is about Ik, i. e., about 20/cbT 1 for typical 
bilayers. However, already at / = 0.5 its magnitude has 
decreased by about 20%. This shows that we have to ex- 
pect a narrowing-down of the hysteresis amplitude com- 
pared to an athermal buckling scenario. 

Upon increasing the adhesion w even further, one will 
reach a critical value wq at which the "back-bending 
branch" of the force-distance-curve touches the vertical 
line / = 0. At this point the tip is being pulled into the 
pore even if there is no force at all. Conversely, neglecting 
barrier complications, this also implies that at the criti- 
cal adhesion energy w — wq an infinitesimal pulling force 
will suffice to unbind tip and membrane, even though the 
adhesion between tip and membrane is greater than zero. 
It is very important to keep this fact in mind if one wants 
to use AFM measurements for the determination of ad- 
hesion energies. 

For w > wq one obtains stable solutions even when 
pulling the tip upwards (where / < 0) The max- 

imum possible force before detachment, /dot, again cor- 
responds to the leftmost point of the back-bend, and it 
increases with increasing w; we will come back to this 
later. Notice that detachment always happens for val- 
ues of ho which are positive, i. e., when the AFM tip 
is still below the substrate level. Contrary to what one 
might have expected, pulling will in this case not draw 
the membrane upwards into a tubular lipid bilayer struc- 
ture (a "tether" ) , which at some specific elongation will 
fall off from the tip and snap back. Rather, the strong 
adhesion pulls the tip far into the pore, and while pulling 
on it indeed lifts it up, unbinding still happens below pore 
rim level. 



2. Strong adhesion energy 

At even larger adhesion energy entirely new station- 
ary solution branches emerge, as Fig. [JJ illustrates for 
w 6 {5, 10, 15} and a = 1. We first recognize the well- 
known hysteretic branch, already seen in Fig. [6] which 
for increasing w extends to much larger negative forces, 
even though the snap-off height ho does only change 
marginally. The shapes of two typical profiles are illus- 
trated in the insets a and b. Notice that this branch 
is always connected to the origin, but for larger values 
of w it starts off into the third quadrant (negative val- 
ues for / and ho)- At first sight it seems that we finally 
get solutions which correspond to a pulled-up membrane; 
however, this region close to the origin corresponds to a 
maximum and is thus unstable. 

a. Overhang branches. Contrary to the hysteretic 
branches, the new branches depicted in Fig. [7J do not 
connect to the origin. This classifies them as a genuinely 



nonlinear phenomenon, since they cannot be obtained as 
a small perturbation around the state / = ho = 0. In 
the first quadrant (/, ho > 0) they all correspond to pro- 
files which show overhangs (see inset c and d). These 
branches had been omitted in Fig. [6l since for weak ad- 
hesion they always correspond to maxima and are thus 
irrelevant. This changes for stronger adhesion, though, 
where they become stable in certain regions (for instance, 
inset c is locally stable). The details by which this hap- 
pens are complicated and will be discussed in more detail 
below. 

Following the new branches to negative forces we see 
that the one for w — 15 loses its overhang around / ps — 6. 
That this can happen continuously is not surprising, since 
within angle-arc-length parametrization there is nothing 
special about the point where \ip\ = 90° (only the shoot- 
ing method might use occurrences of \ip\ > 90° as a po- 
tential termination criterion for integration). 

b. Branch splitting. We also see that (for sufficiently 
large w) there is a point where the hysteresis branch in- 
tersects the new nonlinear branch. There the values of 
/ and ho coincide for both branches, but the detach- 
ment angle a and the total energy of the profile are 
generally different. However, the difference in energy at 
the intersection decreases with increasing w, and around 
w = 15.3 it finally vanishes. At this degenerate point a 
branch splitting occurs, where the connectivity of the two 
branches re-bridges, as illustrated in the lower left inset 
in Fig. [7J Rather than connecting to the origin, the wide 
loop of the original hysteresis branch now joins into the 
overhang branch of the first quadrant, while its bit that 
was connected to the origin now joins into the overhang 
branch in the third quadrant. 

c. Cusps. Figure [5] shows the force-distance curve 
branches for the even larger adhesion energy w = 20. 
This depicts a situation well after the branch-splitting, 
so we recognize the old hysteretic branch connecting with 
overhangs, and the branch connecting to the origin ex- 
tending exclusively in the third quadrant. In contrast to 
Fig. [Jj the line styles in Fig. [8J are chosen to illustrate 
local minima (solid) or maxima (dotted). What imme- 
diately strikes one as surprising is that the profiles at 
/ = —4 belonging to the insets / and h both correspond 
to maxima, even though they sit on both sides of a back- 
bending branch, close to its end (compare this to the 
"usual" scenario at (/ —49, ho » 3.6). Moreover, the 
solution belonging to inset / turns into a local minimum 
for slightly more negative forces, without any noticeable 
features of the branch. How can this happen? 

The explanation is illustrated in the lower left inset 
in Fig. [HI which shows the total energy as a function of 
detachment angle a. Recall that extrema in this plot cor- 
respond to stationary solutions. As can be seen, the en- 
ergy is multivalued, meaning that there exists more than 
one solution at a given detachment angle (these would 
then also differ in their value of their penetration ho). 
But more excitingly, this graph exhibits a boundary ex- 
tremum at a lowest possible nonzero value of a in the 
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FIG. 7: Force-distance curves for a = 1 and w — 5 (solid), 10 (dashed) and 15 (dotted). A thinner line style is used for those 
parts of the curves where the corresponding profiles exhibit overhangs. In the insets a-d profiles for different values of (/,/io) 
are depicted (scaling is h : p — 5 : 3). In the inset on the lower left corner the "branch splitting" is shown as discussed in the 
text (w = 15.0 (dotted line) and 15.5 (solid line)). 



form of a cusp. This is how one can have two successive 
maxima on a curve without an intervening minimum — 
the minimum is simply not differentiable. Hence, there 
is a third solution branch, corresponding to the cusp, at 
which the contact curvature condition from Eqn. (jllbp is 
not satisfied, because this condition is blind to the possi- 
bility of having non-differentiable extrema. Plotting this 
cusp branch also into Fig. [8l we finally understand how 
the switching of a maximum into a minimum happens: it 
occurs at the point of intersection with the cusp branch. 
As the lower left inset in Fig. [5] illustrates, the maximum 
belonging to the solution / joins the cusp-minimum (be- 
longing to solution g) in a boundary flat point, roughly at 
force / = —4.8. For more negative forces this flat point 
turns up, leaving a boundary cusp maximum and a new 
differentiable minimum. Notice that a similar exchange 
happens once more at (/ w —26.4, ho « 1.6). Inciden- 
tally, since at the cusp the contact curvature condition is 
not satisfied, and since this is the only point where the 
adhesion energy w enters, the location and form of the 
cusp branch is independent of the value of w. 

The existence of the cusp branch poses the question, 
whether the solutions corresponding to it are physically 
relevant (at least the ones which are minima). It is not 
so much the lack of differentiability at a cusp minimum 
which causes concern, but rather the fact that it is lo- 



cated at a boundary. Take for instance the E(a) curve 
in the lower left inset of Fig. [8] corresponding to / = —4. 
Now consider a (nonequilibrium) solution which sits on 
the upper branch, somewhere between the solutions g 
and h. To lower the energy, this solution will reduce the 
detachment angle a, thereby approaching the minimum 
at g. But once g has been reached, no further reduction 
in a seems possible, since for smaller values no equilib- 
rium solution exists. The crucial point is that our present 
theory is insufficient to answer what else would be going 
on for smaller a. It could for instance be that there are 
indeed solutions, but they are not time-independent. This 
might be analogous to the well-known situation of a soap 
film spanned in the form of a catenoidal minimal surface 
between two coaxial circular rings of equal radius R. It is 
easy to show that for a ring separation exceeding 1.325 R 
no more stationary solution exists, even though the lim- 
iting profile is in no way singular [58| . However, when 
slowly pulling the two rings beyond this critical separa- 
tion, the soap film does not suddenly rupture. Rather, 
it becomes dynamically unstable and begins to collapse. 
In the case we are studying here, the system drives itself 
to the singular boundary point, and without a truly dy- 
namical treatment it is not possible to conclude whether 
it would remain there or start to dynamically approach 
a different solution. For this reason we do not want to 
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FIG. 8: Force-distance curves for a = 1 and w — 20 (including cusp branch). Solid lines correspond to local minima, dotted lines 
to local maxima. A thinner line style is used for those parts of the curves where the corresponding profiles exhibit overhangs. 
In the insets e - h profiles for different values of (f,ho) are depicted (scaling is h : p = 5 : 3). In the inset on the left lower 
corner the total energy E is plotted as a function of detachment angle a for different forces / (see text for further explanation) . 



overrate the significance of the cusp branch; yet, its exis- 
tence is still important in order to explain the behavior of 
the other "regular" branches, for instance their metamor- 
phosis from maximum-branches into minimum-branches 
or vice versa. 

d. Detachment forces. A measurable quantity in the 
experiment is the detachment force between tip and 
membrane, which is the maximum applicable pulling 
force /(j ct before the tip detaches from the membrane. 
In Fig. [9] this force is plotted as a function of adhe- 
sion energy w for different values of the scaled tension 
a. Starting from a certain threshold adhesion w^^a), 
below which no hysteresis occurs, /dot decreases with in- 
creasing w and exhibits a linear behavior for higher adhe- 
sions. Increasing a also increases the threshold adhesion 
(e. g. u>thr(l/16) = 0.22 compared to Wthr(l) = 1-25). In 
the large- u)-limit fdct/w = finally approaches 

a limit which is independent of k and a and only de- 
pends on the geometry. The elasticity of the membrane 
no longer influences the measurement of the adhesion en- 
ergy - not because the membrane is not deformed, but 
rather because its deformation energy is subdominant to 
adhesion. But for more realistic smaller values of w this 
decoupling does not happen, and adhesion energies can 
only be inferred from the detachment force when a full 
profile calculation is performed. 
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FIG. 9: Scaled detachment force /d c t as a function of scaled 
adhesion energy w for four values of the scaled tension, a — 
yjr, j, 1 and 4 (a increasing from left to right). 



At higher values of a also other qualitative features 
(such as additional instabilities) occur. However, these 
ramifications will not be discussed in the present work. 

e. Tethers. Characteristically, the detachment hap- 
pens at deep indentations (ho close to the maximum in- 
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dentation possible). Long pulled-out membrane tubes 
("tethers"), as they have been studied in the literature 
[29l . f30l . [43l | . are not observed. Even though in our cal- 
culations we find profiles with ho < 0, these solutions 
either correspond to energetic maxima, or they are only 
local minima - with the global minimum at ho > corre- 
sponding to a significantly lower energy. This is a conse- 
quence of the adhesion balance present in our situation: 
upon pulling upwards, it is more favorable for the tip ei- 
ther to be "sucked in" completely or to detach from the 
membrane, rather than forming a long tether. As Fig. [8] 
shows, there is a very small "window of opportunity" at 
/ « — 5 where (locally) stable solutions pulled above the 
surface exist. Yet, their profiles look essentially like the 
ones of inset / or g and show no resemblance to real long 
tethers. Upon increasing the force they become unsta- 
ble, such that the tip either falls of the membrane, or is 
drawn below the membrane plane (notice that there exist 
two minima at / slightly smaller than —5, but both at 
positive indentation) . 

This analysis shows that it appears impossible to pull 
tethers using a probe with a certain binding energy, de- 
spite existing experiments in which tethers of microme- 
ter size were generated [HI, [H, |43[. Consequently, the 
assumption of an adhesion balance does not seem to be 
correct in these cases. Indeed, in these studies the exper- 
imental setup was different (membrane-covered micron- 
sized beads [HI, HU and AFM tips covered by lipid mul- 
tilayers [1611) . In the present situation tethers are also 
observed (26[, but these events are not very reproducible, 
and based on the above calculations we would tentatively 
attribute them to a pinning of the membrane at some ir- 
regularity of the tip. 



V. DISCUSSION 

In the previous sections we have discussed the indenta- 
tion of a pore-spanning bilayer by an AFM tip. We have 
seen that the force-distance curves show a linear behavior 
for small forces in a broad parameter range if the adhe- 
sion between tip and membrane vanishes. Even though 
this is in agreement with recent experiments ([26[, see 
also Fig. O, such a linear behavior is unfortunately too 
featureless to reveal the values of both elastic material 
constants, a and n. 

One way out of this apparent cul-de-sac would be to 
repeat the experiment for different pore radii i? por e while 
keeping all other parameters fixed. Since a and k are the 
same for all pore sizes in that case, it should be possible 
to extract their value from the measured force-distance 
curves. Note that one does not have to fit both parame- 
ters simultaneously if one at first considers a pore where 
the radius is much larger than the characteristic length- 
scale A (see Eqn. |J2])). The corresponding system is in 
the high tension regime and the force-distance relation 
only depends on the surface tension (see Section TlV A[) . 
After determining a from the resulting curve, one may 



subsequently extract the value of k from a measurement 
of a system with smaller pore size. 

The elastic constants can also be obtained by consider- 
ing systems where the adhesion w between tip and mem- 
brane has been increased experimentally. As we have 
seen, the curves change their behavior dramatically for 
ic/0. It should thus be possible to fit two parameters 
to the resulting curves which would yield a local k and 
cr in one fell swoop whereas w can simultaneously be de- 
termined from the snap-on of the tip upon approach to 
the bilayer. The experimentalist, however, has to make 
sure in that case that the line of contact between tip and 
membrane is really due to a force balance as described in 
this paper and not due to other effects such as pinning 
of the membrane to single spots on the tip. In practice, 
this is rather difficult and will be a challenge for future 
experiments. 

One also has to keep in mind that the assumption of 
a perfect parabolic tip is quite simplistic compared to 
the experimental situation. It is probably valid in the 
vicinity of the apex but generally fails further up. Since 
the force-distance behavior close to the depth-saturation 
depends strongly on the actual tip shape, one can only 
use that part of the force-distance curve for data inter- 
pretation where the indentation is still small. To predict 
the whole behavior the exact indenter shape has to be 
known: as long as the situation stays axisymmetric one 
may, in principle, redo the calcutions of this publication 
with the new shape. This is, however, rather tedious and 
therefore inexpedient in practice. 

Our theorical approach does not account for hydrody- 
namic effects although the whole setup is in water and the 
AFM tip is moved with a certain velocity. First measure- 
ments have shown, however, that it is possible to increase 
the velocity of the tip up to 60/xms -1 without altering 
the force-distance curves dramatically [26| . One can un- 
derstand this result with the help of the following simple 
estimate: assume that the tip is a sphere of radius i?ti P 
moving with the velocity v in water. When indenting 
the membrane to a distance d it will also have to over- 
come a dissipative hydrodynamic force i<diss m addition 
to the elastic resistance of the membrane. The energy 
dissipated in this process, -Ediss = -Fdiss^ = 67r7yi? t i p ud, is 
of the order of the thermal energy if typical values are 
inserted (77 = 10~ 3 Pa s, i? t i P = 10nm, v — 60/ims _1 , 
d = 100 nm). This is substantially smaller than the cor- 
responding elastic energy -E e i as t- Complications arising 
from a correct hydrodynamical treatment were thus omit- 
ted here. 

Including adhesion, the velocity of the measurement 
should nevertheless be as slow as possible to ensure that 
the line of contact equilibrates due to the force balance. 
If this is guaranteed, one can also check whether the pre- 
dicted linear behavior between detachment force and ad- 
hesion is actually valid. 
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APPENDIX A: BOUNDARY CONDITIONS 

In this appendix we will explain the origin of the 
boundary conditions j6| and (TTTj) : Eqns. (|6aj) follow sim- 
ply from the requirement of continuity at the pore rim 
and the point where the membrane leaves the tip. Asking 
for a membrane that has no kinks and thus no diverging 
bending energy gives Eqns. (|6bf and (| 1 1 a,|) . 

If the membrane is free to choose its point of detach- 
ment as it is assumed here, an adhesion balance at the tip 
yields another bound ary condition for the contact curva- 
tures ([6c]/(TIb| (see (El Sec. 12, problem 6], [36|, and 
f39l]). In Rcf. [39} a quick derivation can be found for 
the axisymmetric case in the constant height ensemble: 
varying the point of contact changes the energy of the 
free profile but also the energy due to the part at the 
tip. By setting the total variation to zero one obtains 
the well-known contact curvature condition (Eqn. (|6c|) in 
Monge gauge). Observe that this assumes differentiabil- 
ity of the energy as a function of contact point position. 
In the force ensemble an extra term FSh has to be added 
to the variation of the bound membrane. A term that is 
equal and opposite, however, enters the variation of the 
free membrane via the Hamiltonian In total, both 
terms cancel and one again obtains the same condition 
(Eqn. db| in angle-arclength parametrization) . 

The remaining condition (|llc[) stems from the fact that 
the total arclength is not a conserved quantity, which it 
would be if we used a fixed interval of integration. Relax- 
ing this u nphysic al constraint requires the Hamiltonian 
to vanish [3J, 15911 . 



APPENDIX B: CALCULATION OF THE FORCE 
VIA THE STRESS TENSOR 



I) and tangential to the curve. The symbol V_l denotes 
the derivative along I. 

In angle-arclength parametrization, the curvatures are 
given by: K± = —ip, Kh — — sin (i]j)/p, and K — 
—p^/(2p). Eqn. (|B1[) can then be written as 



F = —- 



{[(* 



sin ip 



simp 



-(^-— p)cosV>} • (B2) 



The integrand can be evaluated further by inserting the 
Hamilton equations p0[) and making use of the fact that 
the Hamiltonian ([9]) is zero. One obtains 
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(B3) 



If we now exploit axial symmetry by integrating around 
a circle of radius p = i?i n t, we finally get F = p z + F; the 
momentum p z conjugate to z has to vanish identically 
which implies that the Lagrange multiplier function X z 
is equal to F. This at first maybe surprising result is no 
coincidence at all. In fact, in Ref. [59| it was shown that 
the Lagrange multiplier functions which fix the geomet- 
rical constraints are closely related to the external forces 
via the conservation of stresses. 

Expression (|B1[) can also be translated into "Monge 
gauge" . If we again exploit axial symmetry and integrate 
around a circle of radius p = i?j nt , it reads 



-2irR h 



n(h"{pf h'{pf 
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h"{p) 



p 2 g 



h'(p) 
V9 



h'( P )\> 1 

p^fa' 9 



(B4) 



where g = 1 + h'(p) 2 . Note that the dash denotes deriva- 
tives with respect to p. 

If in particular we choose to evaluate the force at 
flint = -Rporo, the expression (|B4[) simplifies considerably 
to Eqn. ®. 



If the shape of the free membrane is known, the stress 
tensor f a (a e {1,2}) can be evaluated at every point 
of the surface Ef ree . The integral of its flux through an 
arbitrary contour 7 which encloses the tip gives the force 



F 



APPENDIX C: NUMERICAL CALCULATIONS 

The Hamilton equations IjlOp were solved by using 
a shooting method [tioj ]: for a trial contact point c 
Eqns. (|10p were integrated with a fourth-order Runge- 
Kutta method. The value of c determined the contact 
angle a and with it ip, P, Pif>> an d P P at s = s via the 
boundary conditions pip . The integration was stopped 
(Bl) as soon as p was equal or greater than i? p0 re- To reach 



The normal vectors I and n are perpendicular to 7 and 
to each other in every point of the curve. In addition, I is 
tangential to the surface whereas n is normal to it. K± 
and K\\ are the curvatures perpendicular (in direction of 



i? poro exactly one extra integration with the correct step- 
size backwards was performed. Finally, the value(s) of c 
for which ip = at i? p0 rc were identified for fixed param- 
eters F, a, w, etc. 
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If the calculation had been done in the constant height 
ensemble, one would additionally have to check whether 
the correct indentation ho was reached at p = i? p0 re after 
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